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Abstract 

The specific heat of a two dimensional repulsive 
Hubbard model with local interaction is investi- 
gated. We use the two-pole approximation which 
exhibits explicitly important correlations that are 
sources of the pseudogap anomaly. The interplay 
between the specific heat and the pseudogap is the 
main focus of the present work. Our self consistent 
numerical results show that above the occupation 
riT ~ 0.85, the specific heat starts to decrease due 
to the presence of a pseudogap in the density of 
states. We have also observed a two peak struc- 
ture in the specific heat. Such structure is robust 
with respect to the Coulomb interaction U but it is 
significantly affected by the occupation tit- A de- 
tailed study of the two peak structure is carried out 
in terms of the renormalizcd quasi-particlc bands. 
The role of the second nearest neighbor hopping on 
the specific heat behavior and on the pseudogap, is 
extensively discussed. 

1 Introduction 

The one-band Hubbard model with repulsive lo- 
cal interaction [T] is one of the simplest but sig- 
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nificant lattice models which allows to describe a 
rich phenomenology of strong correlated electron 
systems. Such model has become even more impor- 
tant after the discovery of the high temperature su- 
perconductors (HTSC) [2]. Nevertheless, although 
the two dimensional one-band Hubbard model has 
been intensively investigated within different lev- 
els of approximation, there are still open issues, 
for instance, the existence of superconductivity and 
pseudogap [3J 2J [5]. 

The nature of the pseudogap is one important 
key to understanding HTSC. A substantial indica- 
tion of pseudogap can be confirmed through the 
analysis of the specific heat behavior [SI [7] . From 
the theoretical point of view, different approaches 
[SI O [TO] have been used to investigate the specific 
heat of the one-band Hubbard model. However, 
the focus has been on the observed two peak struc- 
ture rather than its connection with the pseudo- 
gap. Monte Carlo simulation and dynamical 
mean field approximation |12| on the two dimen- 
sional one-band Hubbard Model showed evidences 
of a pseudogap. But such calculations were carried 
out only at half filling. 

In this work, we calculate the electronic spe- 
cific heat out of the half-filling (tit < 1, where 
tit = n a + n~ a) for the local repulsive Hubbard 
model. The Green's functions have been obtained 
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where /(w) is the Fermi function and the Green's 
function is 
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in a two-pole approximation [131 114j which takes 
into account essential correlations that allow to 
capture important properties of the pseudogap. In 
addition to the nearest neighbor hopping, the Hub- 
bard model considered in the present work includes 
hopping to second nearest neighbors. Such hopping 
enhances the spin-spin correlations that affect both 
the pseudogap and the specific heat structure. We 
present results for the normal phase, i. e., T > T c 
where T c is the superconducting critical tempera- 
ture. Firstly, the two peak structure of the spe- 
cific heat is analyzed in terms of the renormalized 
quasiparticle band. Then, the pseudogap and its 
interplay with the specific heat, is discussed. 

2 Model and general formula- 
tion 

The model investigated is the two-dimensional one- 
band Hubbard model [T] 

H= ^ tijdlrdjo- + E n i,a n i-a ~ M E 
{(ij))tr ia to 

(1) 

where d irT (di a ) is the fermionic creation (annihila- 
tion) operator at site i with spin a = {f, J,} and 
«i,cr = d\ a di a is the number operator. The quan- 
tity tij represents the hopping between sites i and 
j and ((...}} indicates the sum over the first and 
second-nearest-neighbors of i and \x is the chemical 
potential. U is the repulsive Coulomb potential be- 
tween the d electrons localized at the same site i. 
The bare dispersion relation is 

= 2t[cos(k x a) + cos(fc y a)] + 4i 2 cos^a) cos(k y a) 

(2) 

where t is the first-neighbor and t<i is the second- 
neighbor hopping amplitudes and a is the lattice 
parameter. 

In the two-pole approximation proposed by Roth 
[T31 [13], the Green's function matrix is G (uj) = 
N (cjN — E)~ N in which N and E are the normal- 
ization and the energy matrices, respectively [14] . 

Following the formalism discussed in reference 
[T5] . the energy E = Up- can be written as (N 



where = y/(U e- k + W^f + ^ a U{e % - W % a ). 

In the grand canonical ensemble, the energy is a 
function of the chemical potential E = E(fi(T)), 
where /i changes with the temperature. Therefore, 
the calculation of C(T) defined as C(T) = ff , 
must be performed keeping (n) constant in the 
T — ll plane [5|. Combining equations ([3]), (HJ) and 
C(T) = ff we get C*(T) = F(u)duj where 

F{u) = f'(Lj)g(u) (7) 



in which f'(ui) = ^ 9 qt an< ^ ^ s the Fermi 

function. The function g{uS) is defined as: 

9^) = ^ E E ZiAkW" - <Vg) (8) 

where 

Z ha (k) = (w. ^ + M + e%) u t aU Z^(k)- (9) 

The band shift Wz introduced in equation (|6|), 
is: W kc = nM-n„) 7i T,j<k-q)F a {q), where e(k- 
= E<(i=o)>j#o*o 3 -e i ( fc -^'^ and F a (q) is given 
in terms of (Sj ■ Si), n 0ja = ± £g FuG^e^ 
and m 0ja = ± T^G^ e lk '^ , where F u T{uj) = 
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Figure 1: In the upper panel, the specific heat 
as function of temperature for different values of 
C/.The lower panel shows the chemical potential. 



§ du> f (ui)T (u>) , in which f{to) is the Fermi func- 
tion and T(u) a general Green's function. The 
Green's functions G , l™ m - ) are obtained as in ref- 
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The spin-spin correlation function is 
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where h ih ^ a = a a \fi—{aij- (r nij,-(r+bij,- <r mij- (r )] 
with a a = (1 - b ii(T )/[2(l - bn_ a bn a )\ and 
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3 Numerical results 

The upper panel in figure [1] shows the specific heat 
as a function of temperature for tit = 0.90 and dif- 
ferent values of U . For U = 0.0, the specific heat 
presents a peak in the region of fcsT w 0.55|f|. 
This peak is known as the Schottky anomaly ( 
see ref. [16j and references therein). However, for 
U = 2.0\t\ the Schottky anomaly starts to split giv- 
ing rise to a two peak structure. So, for U = 4.0|£| 
the specific heat shows a peak at low temperature 
ksT ps 0.2|i| and a second peak in ksT ps 1.3|t|. 
For U > 4.0|t|, the specific heat has a two peak 
structure, i.e., the high temperature peak persist 



even in the strongly correlated limit. The lower 
panel in figure [T] shows the chemical potential ii as a 
function of temperature and different Coulomb in- 
teractions. The occupation is tit = 0.90, therefore, 
for U = the chemical potential is found below the 
van Hove singularity. As the temperature increases, 
ix moves to higher energies in order to keep tit un- 
changed. However, when fx reaches the VHS (at 
k B T/\t\ ^ 1.1), the high density of states in VHS 
forces the chemical potential back to lower energies 
to maintain tit unchanged. When U ^ 0.0, a gap 
opens on the density of states near the VHS and \i 
increases with temperature until the Fermi function 
reaches the upper Hubbard band. In this case, for a 
given temperature, \i can be found within the gap. 
If temperature increases more, the upper Hubbard 
band increases its contribution and \i must return 
to lower energies to maintain tit unchanged. The 
behavior of the chemical potential is closely related 
to the specific heat structure, mainly at low tem- 
peratures. We verify that the local minimum in 
C(T) occurs when ll is found within the gap due to 
U. 

We can understand the two peak structure on 
C(T) if we analyze the functions F(u), /'(w) and 
g(oj) defined above. The function g(uj) is shown in 
the left panel in figure[2]for different Coulomb inter- 
actions. The function g{ui) is directly related to the 
renormalized density of states, therefore presents a 
gap due to the Coulomb interaction U. For the 
model parameters considered in figure [21 in gen- 
eral, g(w) is related to the lower Hubbard band 
if ui < 0.50 and to the upper Hubbard band if 
w > 0.50. We verified that at low temperature, the 
function f'(oj) is very concentrated on the chemi- 
cal potential and vanishes quickly for energies that 
are not close to the chemical potential. Moreover, 
for tit < 1, the chemical potential is found to be 
localized in the lower Hubbard band. Then, at low 
temperatures only the lower Hubbard band con- 
tributes to the specific heat and therefore, the low 
temperature peak on the specific heat is associated 
to the lower Hubbard band. On the other hand, at 
high temperatures the function f'(ui) spreads out, 
reaches the upper Hubbard band and the specific 
heat presents a second peak due to the contribution 
of that band. 

Furthermore, as can be verified in figure [TJ the 
position of the second peak in the fc^T axis moves 
to high temperatures as U increases. This occurs 
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Figure 2: In the left panel the behavior of g(uS) for ksT = 0.03|i| and different values of U . The right 
panel shows the function F(u>) for U = 8.0\t\ and different temperatures. In both panels ut = 0.90, 
t = -1.0 and t 2 = 0.3|i|. 




Figure 3: The lines with symbols show the effective spectral weight Zi M (k) and the lines with no 
symbols show the renormalized bands Wi jCT (fc). The upper panels show Z\ M (k) and ui^(k) while the 
lower panels show Z2 jCr (fc) and oj2.a(k). The panels (a) and (b) correspond to U — 8.0\t\ and the panels 
(c) and (d) correspond to U = 12.0|t|. 



because the upper Hubbard band moves to higher 
energies as U increases. Consequently, g(w) asso- 
ciated to the upper Hubbard band moves also to 
higher energies and f'(oS) can reach this region only 
if ksT is large enough. Indeed, we verified that the 
relation between the high temperature peak posi- 
tion and U is fc^T « -j. The right panel in figure 
[5] presents the function F{uj) (defined in equation 
(UJ) for different temperatures. Due to the rela- 
tion between f'(ui) and F(u>), at low temperatures, 
F(u>) is concentrated around the chemical potential 
fi (localized in ^ = 0). Moreover, F(lo) presents a 
negative region near /i that decreases with temper- 
ature resulting in an increasing of the specific heat. 
As a consequence, a peak appears on the specific 
heat at low temperature. For fc^T = 1.00|t|, the 
negative region in F(oj) vanishes and the positive 



region decreases due to the behavior of f'(cd). This 
feature of F(ui) produces a small value for the spe- 
cific heat at fcgT = 1.00|t|. Increasing the temper- 
ature, for instance fc^T = 2.40 \t\, the intensity of 
F(u>) is low for & < 0, and a significant region of 
positive F(uj) emerges for ^ > 5.0. This occurs be- 
cause at high temperatures f'(u>) spreads out and 
reaches the region of g(ui) associated to the upper 
Hubbard band. This behavior of F(u) enhances 
again the specific heat giving rise to the high tem- 
perature peak. For fc^T = 4.0|£|, F(uj) acts as for 
ksT = 2.40 however, it decreases the intensity 
due to the decreasing of the /'(w) intensity. 

Figure [3] shows the renormalized bands u. a ^ and 
the effective spectral weights Zi yCr (k) defined in 
equations ^ and ©• At low temperatures Zi !<7 (k) 
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Figure 4: In the upper panel, the specific heat for 
U = 8.0\t\ and different occupations tit- The lower 
panel shows the behavior of ^ as a function of the 
occupation tit- 



is negative only in a small region next to (§,§)■ 
As temperature increases, this negative region en- 
larges along the direction (§, §) - (0, 7r). However, 
the intensity of Z\ jCr (k) is very small. In^c), the 
negative regions of Zi^(k) at low temperatures are 
greater than in[3|a), but its intensity remains small 
when compared with the intensity of the positive 
regions. 

An important feature observed in lo 1 g is the 
presence of a pseudogap near (0, 7r). This pseu- 
dogap plays an important role relative to the low 
temperature peak observed on the specific heat. In- 
deed, the presence of a pseudogap affects strongly 
Z\, a {k) which is strongly dependent on lo 1 a ^ (see 
equation © ). The opening of the pseudogap in 
(0,7r) moves the renormalized band lo 1 a j: to nega- 
tive energies. As a consequence, Z\_ a (k) is shifted 
to the positive region as can be observed, for in- 
stance, in figure G3c), when ksT = 0.1|t|. In (b) 
and (d), Z2 !<7 (k) is always positive and increases its 
intensity with U. The pseudogap emerges when the 
correlations of antiferromagnetic character associ- 
ated to the correlation function (o, ■ Sj) present in 
the band shift, become strong sufficiently to push 
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Figure 5: The upper panel shows the specific heat 
function of temperature for different values of 
The lower panel shows the spin-spin correla- 
tions (Si ■ Sj) for the same parameters as in the 
upper panel. 

The upper panel in figure U shows the specific 
heat C(T) as a function of temperature for differ- 
ent occupations. Notice that at low occupations 
C{T) is characterized by a peak at low tempera- 
ture. On the other hand, when ut increases a sec- 
ond peak appears at high temperatures. At low tit 
only few electrons reach the upper Hubbard band 
and therefore the high temperature peak on C(T) 
is negligible. The behavior of £ as a function of 
tit is shown in the lower panel of figure [4] This 
results agree qualitatively with those in reference 
[T7] , We verified that for nx > 0.85 a pseudogap 
appears on the anti-nodal points of the renormal- 
ized quasi-particle bands (see the lower panel in 
figure E]) . Such pseudogap suppresses the density 
of states (DOS) on the chemical potential. As the 
specific heat C(T) is directly related to the DOS, 
which in turn depends on the renormalized quasi- 
particle bands, the effects of the pseudogap appear- 
ing also on C(T). In the present work, the function 
F(w) defined in equation Q) associates the renor- 
malized quasi-particle band and the specific heat. 
We verified that the pseudogap suppresses F(u>), 
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~ 0.60. For % > 0.60 the specific heat does 



Figure 6: The upper panel shows the specific heat 
has a function of % and different temperatures. 
The lower panel shows the renormalized quasi- 
particle band <jJ lcr £ for ksT/\t\ = 0.20 and the same 
parameters as in the upper panel. The inset shows 
explicitly the pscudogap Ape- 
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not change significantly. The lower panel in fig- 
ure IB] presents the renormalized quasi-particlc band 
0.20 and three distinct values 
0.30, a pseudogap is observed 

near the antinodal point (0,7r). If increases the 
pseudogap closes as shown for m = 0.50. There 
is a remarkable coincidence between the maximum 
on the specific heat and the maximum pseudogap. 
The analysis of F(ui) shows that it gives the greater 
contribution to C(T) when ^ = 0.30. This occurs 
due to the wide fiat region in w lcr £ which extends 
from (§,7r) to approximately (f , f ) and gives rise 
to a very large F(ui) near the chemical potential. 
Therefore, even that t-i opens a pseudogap which 
suppresses F(yj) on the chemical potential, it also 
produces the flat region (near the chemical poten- 
tial) that overcomes the effect of the pscudogap and 
increases significantly the specific heat. At high 
temperatures the pscudogap closes and also the ef- 
fect of the flat region is suppressed by temperature 
effects. 



4 Conclusions 



below the chemical potential (in uj = 0), resulting 
in a decreasing in C(T) for ut > 0.85. Therefore, 
it can be concluded that the decreasing of the spe- 
cific heat above tit — 0.85 in figure [4] is a clear 
manifestation of the pseudogap. 

The effects of the second-nearest-neighbor ti 
on C(T) are shown in the upper panel in figure 
[5j We observe that the low temperature peak 
on C(T) is more affected by t% than the hight 
temperature peak. The lower panel shows that 
ti enhances significantly the spin-spin correlations 
(Si ■ Sj) which modifies the renormalized quasi- 
particlc bands and consequently the function F(uj) 
and therefore, C(T). At low temperature \ (Si ■ Sj}\ 
is stronger and then its effects on C(T) are more 
evidenced than at high temperatures. This is the 
main reason why the low temperature peak on 
C(T) is more intensively affected by t^. 

The upper panel in figure [5] shows the specific 
heat as a function of the second nearest neighbor 
hopping amplitude & at different temperatures. 
For ksT/\t\ = 0.20 the specific heat presents a 
maximum for ^ 0.30 and then it decreases until 



In summary, the analysis of the two peak structure 
of C (T) in terms of the renormalized quasi-particle 
bands allowed us to investigate also a very impor- 
tant feature present in the region of Ut ^ 0.85, i. 
c. the pscudogap region. We observed that above 
riT ~ 0.85, the specific heat decreases signaling the 
pscudogap presence. In reference the authors 
suggest the presence of a pseudogap at half filling of 
the Hubbard model with Monte Carlo simulation. 
Here, we confirm the pseudogap existence also out 
of the half filling. 
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